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WMAP 3yr data with the CCA: anomalous emission and impact of 
component separation on the CMB power spectrum 
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ABSTRACT 

The Correlated Component Analysis (CCA) allows us to estimate how the different diffuse 
emissions mix in CMB experiments, taking also into account complementary information 
from other surveys. It is especially useful to deal with possible additional components for 
which little or no prior information exists. An application of CCA to WMAP maps assum- 
ing that only the canonical Galactic emissions (synchrotron, free-free, and thermal dust) are 
present, highlights the widespread presence of a spectrally flat "synchrotron" component, 
largely uncorrected with the synchrotron template, suggesting that an additional foreground 
is indeed required. We have tested various spectral shapes for such component, namely a 
power law as expected if it is flat synchrotron, and two spectral shapes that may fit the spin- 
ning dust emission: a parabola in the log S'-log v plane, and a grey body. If the spatial dis- 
tribution of the additional ("anomalous") component is not constrained a priori it is found to 
be always tightly correlated with thermal dust, but the correlation is not perfect. Quality tests 
applied to the reconstructed CMB maps clearly disfavour two of the models. The CMB power 
spectra, estimated from CMB maps reconstructed exploiting the three surviving foreground 
models, are generally consistent with the WMAP ones, although at least one of them gives a 
significantly higher quadrupole moment than found by the WMAP team. Taking foreground 
modeling uncertainties into account, we find that the mean quadrupole amplitude for the three 
"good" models is less than Icr below the expectation from the standard ACDM model. Also 
the other reported deviations from model predictions are found not to be statistically signif- 
icant, except for the excess power at I ~ 40. We confirm the evidence for a marked north 
south asymmetry in the large scale {I < 20) CMB anisotropics, which is stable with respect 
to the foreground parametrization we adopted. We also present a first, albeit prehminary, all 
sky map of the "anomalous emission". 



1 INTRODUCTION 



The Wilkin s on Microwave Anisotrop y Probe (WMAP, 
iBeiinett et al] l2003l : ISpergel et all l2006f) data are having a 
dominant role in determining the reference cosmological model. 
It is therefore important to get an in depth understanding of all 
the uncertainties in the data analysis and of their effects on the 
estimates of the angular power spectrum of the Cosmic Microwave 
Background (CMB), from which the values of the fundamental 
cosmological parameters are derived. 

In this paper we present a new investigation, exploiting a novel 
technique, of the effect of foreground contamination. CMB power 
spectrum analyses are normally done removing a substantial frac- 
tion of the sky (~ 25%) around the Galactic plane where Galac- 
tic emissions are stronger. However, the uncertainties in the under- 
lying foreground emission make it possible that a relevant contri- 
butio n from them remains, even when large sky cuts are consid- 
ered, [sios^^ Seliak|j2004h developed a formalism to incorporate 
foreground uncertainties and the effect of sky cuts in the estimates 
of errors on multipole moments. ISlosar et al. (200 4) marginalized 
over the foreground templates produ ced b y iBermett et alJ 1 2003 ) 
from WMAP data. lO'Dwveil ( l2005h and lO'DwveretZl 120041) 
have shown that it is possible to include the foreground compo- 
nents in a self-consistent fashion within the statistical framework of 



their approach to power spectrum estimation, based on a Bayesian 
statistical analysis of the data utilizing Gibbs Sampling. 

Here we deal with foregrounds applying to the three- 
year WMAP data the Correlated Co mponent Analysis (CCA, 
iBedini et al1l2005l ; iBonaldi et al.ll2006l ), which exploits a second 
order statistics of the data to estimate the spectral behaviour of 
the mixed components directly from the WMAP data. The method 
also allows us to incor porate the add it ional i nformation from other 
surveys, such as the iHaslam et alj (Il982h 408 MHz map, the 
IRAS/COBE dust emission map ("Finkbeiner et al."l999'), and the 
extinction corrected Ha map (Dickinson et al. 2003), as a tracer of 
free-free emission. At variance with most alternative techniques, 
CCA allows us to perform an essentially all-sky component sepa- 
ration: only a narrow strip around the Galactic equator (\b\ < 3°, 
~ 5% of the sky) needs to be removed dBonaldi et alj200q) . 

Of course, we need to specify a priori which foreground com- 
ponents are contributing to the maps and describe their frequency 
spectra with a small number of parameters, as specified is § 2, and 
this brings in the issue of the possible presence of one more Galac- 
tic component, in addition to the well established synchrotron, free- 
free, and thermal dust emissions. In § 3 we describe the estima- 
tion with CCA of the mixing matrix, exploited to reconstruct the 
maps of astrophysical components using an inversion algorithm 
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that, in our case, is tiie harmonic Wiener Filter, as described in § 4. 
The properties of the reconstructed foreground components are dis- 
cussed in § 5, and those of the CMB in § 6 where we analyze the un- 
certainties on the derived CMB power spectrum ensuing from those 
on foreground properties, especially in relation to the reported devi- 
ations from predicti ons of the "concordance" cosmological model 
dSpergel et alj2006l) . The main conclusions are summarized in 5 7. 



2 DIFFUSE COMPONENTS IN THE WMAP BANDS 

We summarize here the frequency spectra of diffuse components 
in terms of antenna temperature. The CMB blackbody spectrum is 

{hv/kTcMBY ex.p{hu/kTcMB) 



Table 1. WMAP channels 



7a,cmb('^) ~ 



{exp{hu/kTcMB) - 1)2 



(1) 



where h and k are the Planck and the Boltzmann constant, respec- 
tively, andTcMB = 2.726 K. 

The free-free emission spectrum is also quite uniform and well 
approximated, in the frequency range of interest, by: 



(2) 



It is difficult to obtain a spatial template of this emission since syn- 
chrotron dominates radio surveys and thermal dust dominates far- 
IR surveys. A good tracer of free-free emission are Ha maps, but 
they need to be corrected for dust extinction, and this correction 
adds a significant uncertainty (Dic kinson et al.ir2003) . 
The synchrotron spectrum is usually modeled as: 

T'A,synch(j^) oc (3) 

Its spectral index, Ps, is known to vary across the sky, re- 
flecti ng the variatio ns of the energy index of rel ativistic elec- 
trons ( Rvbicki & Lig htman 1 97^ . The WMAP team teennett et al] 
l2003l : iHin shaw et al. 2006) find steeper values (/3s ~ 3) at high 
Galactic latitudes and flatter values (/3s — 2.5) close to star- 
forming regions, especially in the Galactic plane. The mean value 
is 2.65 between 408 MHz and the K-band, and 2.69 between 408 
MHz and the Q-band. A spectral steepening is expected at high mi- 
crowave f requencies, as a consequence of electron ageing effects 
( iBandav & Wolfendale 1991). The Haslam et al. ( 1982) 408 MHz 
map is commonly used as a synchrotron template. 

The spectrum of the thermal (vibrational) emission from dust 
grains in the interstellar medium is conveniently described by a 
grey body: 



7A,dust(l^) OC 



,,/3d + l 



exp(/ii//fcrdus 



(4) 



with Tdust ^ 18 K and /3d ~ 1.67 at mm wavelengths 
( lFinkbeineretal.ll 19991) . Both Tdust and (3^ are expected to vary 
across the sky due to variations of the interstellar radiation field 
and of the dust grains composition. Measurements of /3d generally 
lie in the range 1.5 ^ /3d ^ 2.0, but an even broader range has 
been observed by the PRONAOS experiment (Dupac et al. 2003): 
0.8 < /3d ^ 2 4. The WMAP data only weakly constrain the dust 
spectral index dffinshaw et al.l2006l) . The most frequently used spa- 
tial template of therma l dust emission is the one worked out by 
iFinkbeiner etai] d 19991) combining IRAS and COBE maps. 

In recent years evidence has been reported for an additional, 
thermal dust correlated component called "anomalous microwave 
emission", that may dominate, at least in some Galactic regions, 
in the 20-40 GHz range, where it has a spectrum similar to syn- 
chrotron. It may be originated by tiny, fast spinning, dust grains 
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<0.23 


o-o(mK) 
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1.45 
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6.50 


mean pixel rms (mK) 
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0.07 


0.07 


0.09 


0.1 




dDraine & Lazarianlll95 


8l:lde Oliveira-Costa et al.ll 19991 l2002h 



dications of this emission were found by analyses of CMB ex 



periments like S askatoon Jde Oliveira-Costa et alj "l997h. 19 GHz 
(de Oliveira-Cos ta et al.l 1998h. Ten erife Oliveira-Costa et al.1 



1l999. ,2002i : iMu kheriee et al. 2001); Python V dMukheriee etai] 



2003). Although Bennett et al. (2003) concluded that spinning dust 
emission contributes less than 5% of the WMAP Ka-band an- 
tenna temperature, evidence of the anomalous emission were un- 
covered combining WMAP data with other measurements, espe- 

^ 1 ' ' 1 

cially at lower frequencies (Lagache 2003; de Oliveira-Costa et alj 
[2004; Finkbeiner 2004; Watson et al. 2005; Davies et al. 2006|). 
Moreover, the polarization properties of the dust-correlated low- 
frequency component differ from those of the component well cor- 
related with the Haslam template, suggesting a different emission 
mechanism dPage et aT 2006). Nevertheless, the exact nature, the 
spectral properties, and the spatial distribution of this foreground 
remain uncertain. There is not general agreement even on its exis- 
tence; a combination of free-free emission and of strongl y self ab- 
sorbe d synchrotron could also account for the data ( Hinsh aw et al.l 
l2006l) . 

Point sources are another important foreground component, 
but their contribution to the power spectrum is important only on 
scales ^1°, while in this analysis we are mostly interested on 
larger scales, where diffuse emissions dominate. For the present 
analysis we will simply need to mask the brightest sources, as de- 
scribed in the following sections. 



3 MIXING MATRIX ESTIMATION WITH CCA 
3.1 Outline of the method 

As usual, we express the data vector x in each pixel as: 

X = Hs + n (5) 

where H is a M x mixing matrix, s is the A'^-vector of sources 
(components) and n the M-vector of instrumental noise (Al is then 
the number of independent maps used in the analysis). The generic 
element h^c of the mixing matrix is proportional to the spectrum 
of the c-th source at an effective frequency within the d-th sensor 
passband. Using eq. dU we implicitly assume that the spatial pat- 
tern of the physical processes is independent of frequency and that 
the effects of the telescope beam has been equalized in all channels. 

Given a generic signal X, defined in a two dimensional space 
with coordinates rj), the covariance matrix of this signal is: 



Cx{r, ^) = {[X(r, i^) - [X(C + r,r^ + ^)- fif) 



(6) 



where (...) denotes expectation under the appropriate joint proba- 
bility, /I is the mean vector and the superscript T means transpo- 
sition. Every covariance matrix is characterized by the shift pair 
(r, ?/;), where r and tp are increments in the ^ and 77 coordinates. 
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From eq. ([5J we can easily derive a relation between the data 
covariance matrix C^, at a certain lag, the source covariance matrix 
Ca at the same lag, the mixing matrix H, and the noise covariance 
matrix Cn' 



C,(r, V) = HC,(r, V')H' + C„(r, V). 

The covariance matrix can be estimated from the data as: 



(7) 



CAr,^) = ^^[x(C,r,)-/.,][x(C + r,r, + 0)-M.]^ (8) 

where Np is the number of pixels sampling the data. Given a noise 
process, we can model the noise correlation matrix C„: for exam- 
ple, if noise can be assumed signal-independent, white and zero- 
mean, for (r, i/)) — (0, 0) C„ is a diagonal matrix whose ele- 
ments are the noise variances in the measured channels, while for 
(r, tp) 7^ (0, 0) Cn is the null AI x M matrix. If the noise process 
deviates significantly from this ideal model, C„ can be computed 
using Monte Carlo on noise maps in the same way we did for C^. 

Once Cx and C„ are known, eq. Q can be used to identify 
the mixing operator H. The strategy of CCA is to parameterize the 
mixing matrix to reduce the number of unknowns and to take into 
account enough nonzero shift pairs (r, ip) to estimate both H and 
Cs- To solve the identification problem we perform the minimiza- 
tion: 



(r(:,:)) =argmin^ |1 H(r)a[E(T, ^)]H^(r) 



(9) 



T,1p 



-C4r,7/>)-C„(r,V) II 



where F is the vector of all parameters defining H and :) is 
the vector containing all the unknown elements of matrices Cs for 
every shift pair. 

The main product of CCA is an estimate of the mixing ma- 
trix H. The result is defined apart from the normalization of the 
different signals at a given frequency. Therefore we estimate a nor- 
malized mixing matrix, obtained by assigning value 1 to all the 
elements of an arbitrarily chosen row. This matrix can be used to 
perform the source reconstruction with standard inversion methods. 
In our case, this is done by harmonic Wiener Filtering, as we will 
describe in §|4l 



3.2 Input data 

The basic data are the three-year WMAP maps in the K, Ka, Q, V 
and W bands, whose main characteristics are summarized in Table 
[T] The values of ctq reported for every channel allow us to calcu- 
late a noise map at each frequency as a = lyoN^^J , where TVobs 
is the map of the effective number of observations, provided with 
the data. The mean rms per pixel relative to each channel at full 
resolution is also reported. 

The maps were then smoothed to the common resolution of 
1° (Bennett et al. 2003). Since we do not have an exact formula 
to associate an rms level to each map after smoothing, we therefore 
simulated a set of ten WMAP noise maps at each WMAP frequency 
using the corresponding cro and iVoba- We then smoothed them to 
the 1° resolution, and computed the mean rms for each channel. 

We assumed the noise process to be gaussian with variance 
equal to the mean of the all-sky variances of the simulated noise 
maps. Then C„{T,tp) is a diagonal matrix whose elements are 
these empirical variances for {T,tp) — (0,0), and the null ma- 
trix for (r, ?/;) 7^ (0,0). In eq. ^ we then treat C„(r, ?/>) as a 



diagonal matrix whose elements are these empirical variances for 
(r, ^p) — (0, 0), and as null matrices for (r, ip) 7^ (0, 0). 

Strong point sources can contribute substantially to the mean 
surface brightness in some pixels and, since their frequency spec- 
tra are generally different from those of diffuse components, they 
may bias the mixing matrix estimates. Source fluxes are highly di- 
luted in the relatively low resolution maps used here and we find 
that the point source effects become negligible once we mask a re- 
gion of 1° radius around each source contained in the three-year 
point source catalogue provid ed with the WMAP data. In addi- 
tion, based on the analysis bvl BonaldietalJ 1 I2OO6I) , we applied a 
Galactic cut of ±3°. Finally, we masked few highly contaminated 
regions, namely CenA, the Large Magellanic Cloud, p Oph, Orion 
A, Orion B and Tau A. 

To get a better leverage for foreground differentiation we com- 
plemented the WMAP data with a "thermal dust" and a "syn- 
chrotron" map. The former was obtained extrapolating to 850 GHz 
the dust map bv lSchleg el et al. using the best fit model of 

iFinkbeiner et alj 1 19991) . The latter is based on the iHaslam et al.l 
a 408 MHz map. Even if the Haslam map is often assumed 
to be pure synchrotron, contributions from the free-free emission 
can b e important, especially on the Galactic plane (Paladi ni et al.l 
l2005h . Our "synchrotron" map was then obtained by subtracting 
from the Haslam map the free-free cont ribution estimated from the 
Ha map corrected for dust absorption JPickinson et al.ll2003h . 



3.3 Analysis 

To reduce the number of unknowns to a manageable level we need 
to parameterize the mixing matrix H. This means we have to ex- 
ploit our knowledge of foregrounds to model the WMAP data in 
terms of a small number of parameters. On the other hand, as men- 
tioned in § 1, we are not even sure about the number of foreground 
components that need to be taken into account. 

Also, foreground spectral parameters vary across the sky. 
Therefore we would ideally want a map of the values of such pa- 
rameters with sufficiently high resolution. However, CCA needs 
a large n umber of inde p enden t pixels to derive them. As dis- 
cussed by iBonaldi et al] ( l2006h a good compromise is to apply 
CCA to patches of about 1500 deg'^. We worked with patches 
of {Al,Ab) = (50°, 30°) on the Galactic plane, and increased 
Al with increasing |6| as necessary to roughly preserve the 
pixel number (and the patch area). We centered our patches 
at longitudes Ic = {0°, 40°, 80°, ...320°} and latitudes be = 
{0°, ±15°, ±25°, ±35°}. 

The different sky patches intersect each other, but the overlap 
is small enough for them to be taken as independent of each other. 
In fact, the covariance matrices [eq. l[8]l] of adjacent patches are all 
substantially different from each other. 

For any given input model CCA has then yielded the mixing 
matrix for each of the 63 patches. As discussed in the following 
subsection, the first use of these results was to give us some guid- 
ance in the selection of input models. In this analysis and through- 
out the whole pap er we used the HEALPix tessellation scheme 
jGorski et alfcoOSh Fl 



see http://healpix.jpl.nasa.gov I 
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Figure 1. Distribution of the synchrotron spectral index (is obtained by 
CCA for model Ml. 



Figure 2. Map of the synchrotron spectral indices recovered by the CCA 
for model Ml. 



33.1 The standard foreground model, Ml 

We started our analysis adopting the most conservative model 
(which will be referred to as Ml), containing the standard mixture 
of CMB, synchrotron, thermal dust and free-free. This model has 
only one free parameter, the synchrotron spectral index 13s . In fact, 
the spectral behaviour of the CMB and of the free-free emissions 
are known [eqs. (T} and l[2j], and the thermal dust spectrum is only 
very poorly constrained by WMAP data, so that the results are very 
weakly dependent on the assumed dust temperature and emissivity 
index. We adopted the commonly used values Tdust = 18 K and 
/3d = 1.67 [eq. (0]. 

The distribution of f3s obtained with CCA, shown in Fig. [T] 
presents two peaks. To determine position and width of the peaks 
we divided the distribution in two, /3a ^ 2.4 and /3s > 2.4, and 
fitted them with gaussian distributions. One peak is at /3s ~ 2.7 
and has a dispersion a ~ 0.2, which is roughly what is expected 
for synchrotron. The second one, at /3s = 2.3852, is extremely nar- 
row (T ~ 0.0004, hinting at a different component. We explicitly 
checked that this distribution is unaffected by different choices of 
the thermal dust parameters, within the observed ranges (§ 2). This 
immediately begs the question, what is this second flatter compo- 
nent. 

In Fig.|2]we show the map of the recovered synchrotron spec- 
tral indices. Mean values are adopted in the regions where different 
patches overlap. The flatfish component is mainly located at low 
latitudes, more than ~ 40° away from the Galactic center, and is 
not well correlated with the synchrotron template. In the next sub- 
section we investigated different possible spectral shapes for this 
"anomalous" component. 



3.4 Models for the "anomalous" component 

As pointed out by many authors (e.g. de Oliveira-Costa et al. 
2004; Watson et al. 2005), the spectral behaviour of the anomalous 
emission clearly differs from those of free-free and synchrotron 
at 10-15 GHz (where we don't have all-sky maps), but is sim- 
ilar to them in the 20-40 GHz range. This means that any at- 
tempt to estimate simultaneously the spectral parameters of syn- 
chrotron and of the anomalous emission from WMAP data is li- 
able to strong aliasing effects. We have therefore chosen to esti- 
mate the spectral parameters of the anomalous emission keeping /3s 
fixed, and repeating the analysis for several values of it in the range 
2.8 ^ /?s ^ 3.1. Specifically we ran CCA over all sky patches for 
/3s = {2.80, 2.86, 2.92, 2.98, 2.04, 3.10}. 



Anomalous emission models 




Figure 3. Frequency spectra of the anomalous emission for the models M2, 
M3 and MS. The curves are normalized at the K band. For each model, the 
shaded area shows the range of values allowed for the esitmation, while the 
symbols refer to the central values. 



As for the anomalous emission spectrum, we have investigated 
the models summarized in Table|2]with their free parameters, and 
described in the following. The proposed frequency scaling for the 
anomalous emission are also shown in Fig. [3] 



3.4.1 Model M2 



Th e results b y'de Oli veira-Costa et al] ( l2004h . lWatson et al.l ( l2005h . 
and lOavies^ et al. ( 200^ suggest that the spectrum of the anomalous 
emission for v > 20 GHz may be represented by a parabola in the 
(log u, log 5) plane with i^max = 20 GHz: 



logrA,x(j^) = const + 



-meo log Un 



log(i^max/60 GHz) 



21og(f„,ax/60GHz) 



- 2 log 1/ + 



(log 



(10) 



with 1/ in GHz. The free parameter, meo, is the angular coefficient 
at the frequency of 60 GHz in the (log log S) plane. 

Running CCA for the above set of values for /3s we found 
that the mean values of meo over the sky patches are in the range 
3.8 ^ meo ^ 4.5, and correlate with /3s. The linear best fit rela- 
tion, shown in Fig.|4] is: 



meo 



(2.1101 ± 0.0005)/3s - (2.073 ± 0.002). 



(11) 
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Table 2. Summary of the models 



Model 


investigated component 




parametrization 


free parameters and ranges 


Ml 


Synchrotron 




eq.(D 


2.0 ^ /3s < 3.5 


M2 


Anomalous emission 


eq. 


(IO);fmax = 20 GHz 


1.0 ^ meo s; 5.0 


M3 


Anomalous emission 




Ta.xH <x 


2.0 ^ ft ^ 2.6 


M4 


Anomalous emission 




^l]/[exp(/ii//A:T:,) - 1]; = 2.4 


0.2K ^ Tj; ^ 0.7K 


MS 


Thermal + Anomalous dust 


eq. 


{Il;t'max = 20 GHz 


1.0 < meo ^ 5.0, 0.1 < Reo ^ 2.0 



The corresponding s pectral shape is comp atible with the anomalous 
emission detected bv lPavies et alj | |2006|) . In particular, our scaling 
between the K and Ka bands is almost the same they find, and the 
one between Ka and Q is very similar. Our model falls down more 
rapidly then theirs at higher frequencies, but is still inside their error 
bars. 



3.4.2 Model M3 

Alternatively, the anomalous emission might be interpreted as flat- 
spectrum synchrotron, possibly highly self absorbed, associated 
to strong magnetic fields local to star forming regions (and thus 
dust-correlated). To investigate this possibility we parameterize this 
component as a power law, with a spectral index /32 • 

As before, we obtained the distribution of P2 over the sky 
patches for each value of /3s. Such distribution turned out to be 
quite narrow (dispersion a = 0.003) and independent of Ps. The 
mean value of the spectral index 



(/Ja) = 2.144 



(12) 



turns out to be remarkably close to that of free-free, hinting at the 
possibility of aliasing effects. We will come back to this possibility 
in §111 



3.4.3 Model M4 

In this c ase we adopt for the a nomalous component the spectrum 
used by iTegmark et alj j2000l) : it is a grey body, of the form of 
eq. with temperature Tx — 0.25 K and emissivity index /3x = 
2.4. We checked that a simultaneous estimation of both f3x and Tx 
was not feasible, as the algorithm failed to converge. Then, we fixed 
Px = 2.4 and allowed Tx to vary in the range 0.2 K^Tx ^0.7K. 

As it can be seen looking at Fig. [3] M2 and M4 yield quite 
different frequency scalings. Infact, the latter has been proposed to 
fit typical Draine & Lazarian ( 1998) spinning dust models, but it 
is not suited to reproduce the results bv lPavies et al. Given 
the slope between the K and Ka bands, M4 goes down more rapidly. 

The results (Fig.O can be summarized as: 



Tx = 



0.43 ± 0.01 K if/3s^3.0 
0.47 ± 0.01 K if/3s>3.0. 



(13) 



The spectrum rises from K to Ka band, and slowly decreases going 
to higher frequencies. The slope between 60 and 94 GHz is consis- 
tent with what found for M2. 



3.4.4 Model MS 

We tested the possibility that there is a perfect correlation between 
thermal dust and anomalous emission, so that they can be treated 
as a single component parameterized as: 




2.70 2.80 2.90 3.00 3.10 3.20 

Figure 4. Mean estimated values of mgo, with their errors, vs assumed 
values of the synchroti'on spectral index /3s (model M2). 




0.40 

2.70 2.80 2.90 3.00 3.10 3.20 



Figure 5. Mean estimated values of Tx , with their eiTors, vs assumed values 
of the synchrotron spectral index /3s (model M4). 



7A,dust+x(i^) OC TA,dust(j') + Rao ■ 2a,x(j'), 



(14) 



where TA.dust is given by eq. l[4j and Ta,x by eq. l llOt . The factor 
Rao quantifies the relative intensity of the two emissions at a fre- 
quency of 60 GHz. The free parameters of this model are Reo and 
mao, whereas the other ones are fixed as in M2. 

The recovered values of the parameters are shown in Fig. [6] 
The values of mao are within the range found for model M2, which 
did not assume any correlation with thermal dust, but the correla- 
tion with Ps is now weaker. The linear best fit relations of meo and 
-Reo with /3s (solid lines in Fig.|6) are: 



Reo = (-1.4 ± 0.4)/3s + (4.4 ± 1.2) 
meo = (1.3 ± 0.3)/3s + (0.3 ± 0.1). 



(15) 
(16) 
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2.70 2.80 2.90 3.00 3.10 3.20 




2.70 2.80 2.90 3.00 3.10 3.20 



Figure 6. Mean estimated values of Rao and meo, with tlieir errors, vs 
assumed values of the synchrotron spectral index Ps (model MS). 



This is again consistent with lPavies et alj |200(^. both for the slope 
meo and the intensity ratio of the anomalous component respect to 
thermal dust i?60- 



4 SOURCE RECONSTRUCTION WITH THE WF 

Once the mixing matrices of the diffuse components in the WMAP 
data were estimated for the different models, we performed the 
source reconstruction with the harmonic Wiener Filtering. In the 
harmonic space the problem stated in eq. (O simply becomes: 



(17) 



where the vectors Xi,„, s;^, and ni,„ contain the harmonic coef- 
ficients of channels, sources and instrumental noise, respectively. 
Using a linear approach to component separation, an estimate of 
the vector sjm, say sim, can be obtained as: 



Sim H 



(18) 



where the reconstruction matrix W|^' depends on the estimate H 
of the mixing matrix. In the case of Wiener Filtering, we have: 



W 



(0 



(19) 



where Cs and Cn are the source and the noise power spectra, 
respectively. If the source and the noise processes are uncorre- 
lated, Cs and Cn are diagonal matrices. Even in this approxi- 
mation, however, reconstructing the components with WF requires 
the knowledge of the diagonal elements Ci''(i) = Cs'''(i,i) and 



Figure 7. Apodized mask used to preprocess the input maps for WF. The 
minimum (dark grey) and maximum (light grey) values are and 1 respec- 
tively. 



Cn\i) = Cn^'^(i, i). The latter can be easily modeled; in the case 
of uniform noise, we have: 



(20) 



where Ni is the number of pixels, and Oi is the rms pixel noise. On 
the other hand, the power spectra of the components, Ci'\ are not, 
or on l y poorly, known a prio ri. A common approach tHobson et al.l 
Il998l : Istolvarov et al.ll20o3) is to start from a rough estimate of 
the source power spectra ci'' (i) and get the reconstruction matrix 
through eq. l |19l l. The latter allows us to get the estimated com- 
ponents through eq. l |18l l. After the first reconstruction, we obtain 
Cf\i) as the unbiased estimators of the power spectra of the re- 
constructed sources and we update w|^' . The procedure is iterated 
until convergence on the power spectra is reached. 

To derive the initial power spectra, we performed a first, low- 
resolution, reconstruction through eq. J18t exploiting the recon- 
struction matrix: 



W 



(0 



(H^(C, 



''H)~'H^(C, 



(21) 



which corresponds to a maximum likelihood analysis. Note that in 
this case, besides the estimate of the mixing matrix H, recovered by 
CCA, only the prior knowledge of the noise power spectra Cn*^'^ 
is required. We then use as initial for the WF the power 

spectra of these reconstru cted sources ext rapolated to small scales 
assuming behaviour (' ZavariselbOOTh . After convergence has 
been reached, dust, free-free and, if present, anomalous emission 
behave as between multipoles ~ 10 and few hundreds, while 
synchrotron behaves as 

The approach just described uses a minimum number of pri- 
ors: in practice, the components are identified only by their fre- 
quency scalings. The lack of prior covariance information certainly 
limits the quality of the reconstruction, but in this way we avoid 
biasing the results. 

4.1 The WF data set 

For the matrix inversion we used the WMAP K, Ka, Q, V, and W 
bands at their original resolution. Since we now work in the har- 
monic domain, a frequency-dependent beamwidth can be naturally 
accounted for. We assumed circularly symmetric gaussian beams 
with the FWHM values reported in Table[T] Even if this may not be 
a good approximation of the real beams, the effect of beam asym- 
metries is unlikely to be relevant on the relatively large angular 
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scales we are considering here. The noise processes have been as- 
sumed to be uniform over the whole sky: to compute the noise 
levels we averaged the standard deviations of a set of simulated 
WMAP noise maps. 

To take advantage of all the available information we comple- 
mented the WMAP data with the synchrotron, dust and free-free 
maps adopted for the mixing matrix estimation (see §[3]l. These 
maps have resolutions of 60', 5' and 60', respectively. We have al- 
lowed for the effect of variations of spectral indices of synchrotron 
and dust across the sky (while the inversion is done adopting ev- 
erywhere the mean values obtained with CCA) by attributing to 
the corresponding maps uncertainties much larger than the mea- 
surements errors. In particular, we computed the standard devia- 
tion of each prior map outside the apodized mask of Fig. |7] and 
assumed a constant rms level of 20% for the 408 MHz synchrotron 
map and of 10% for the 850 GHz thermal dust map. The spread of 
spectral indices may induce even larger uncertainties, but we have 
checked that increasing the uncertainties by up to a factor of two 
has a negligible impact on the reconstruction of the CMB and a 
small effect on that of foregrounds. This problem does not affect 
the free-free for which we have adopted the nominal rms error of 
7% (Dickinson et al. 2003). 

Algorithms working in the spherical harmonics space require 
data maps defined over the whole sky. On the other hand, cutting 
out highly contaminated regions is generally necessary to recover 
the CMB map. But the spherical harmonic basis is no longer or- 
thonormal over the cut sphere and correlations are introduced be- 
tw een harmonic coefficie nts at different (I, m) modes. As shown 
by IStolvarov et alj ( l2005h . one way to overcome these problems 
is to consider the cut areas as an extreme case of anisotropic noise, 
where the noise becomes formally infinite. In this way we can work 
with all sky maps but the solution is not constrained in the omitted 
areas. Similar results are obtained, with less computational prob- 
lems, by filling the masked areas with a noise realization at the 
nominal level for the considered channel (Stolyarov, private com- 
munication). To reduce edge effects the maps were weighted with 
the apodized mask shown in Fig. [T] This mask excludes the strip 
at |6| < 3°, the LMC, the Cen A region, and regions of 1° radius 
centered on the brightest point sources. Obviously, all pixels which 
have been modified will be excluded from the analysis of the output 
maps. 

4.2 Analysis 

The CCA gives, for each model, the distribution of spectral param- 
eters found in different sky patches. These distribution reflect (be- 
sides estimation errors) the spatial variability of the spectra. On the 
other hand, the harmonic Wiener Filtering is applied using a sin- 
gle mixing matrix over the full sky, i.e. assuming spatially uniform 
spectra. To investigate the effect of this simplification we generated 
a set of 30 mixing matrices for each model, drawing the spectral 
parameters at random from CCA distributions. For each matrix, we 
performed a full-sky source reconstruction up to / = 300. 

In the cases of models M2 and M5, we adopted, for 13s, the 
same distribution of Ml, and derived the value of meo or of Rqo 
and mgo from eq. dl It or from eqs. dlSb and ( 116b . respectively. In 
the case of M3, we used the same uniform distribution for Ps and 
a Gaussian distribution for P2 ■ For M4 we used the same uniform 
distribution for (is and derived Tx according to eq JI3t . 

A particular case is given by the standard foreground model. 
Ml. The results obtained with CCA for this model indicate that, if 
no additional component is present, we cannot assume a uniform 



synchrotron spectral index all over the sky. Thus, generating a set 
of synchrotron spectral indices according to the distribution shown 
in Fig. [T] would not be meaningful. For this reason we used this 
model for reference, adopting for the synchrotron spectral index a 
uniform distribution in the range 2.8 ^ j3s ^ 3.1. 



5 RECONSTRUCTED COMPONENTS 

In Fig. [8] we show the rms in the WMAP bands of the components 
reconstructed using each model. Only pixels outside the apodized 
mask of Fig.Qhave been taken into account. 

On this plot, the CMB component shows tiny variations from 
one iteration to another (the shaded regions are too narrow to be vis- 
ible) and also from a model to another, suggesting that the choice 
of the foreground model, among those considered here, has little 
influence on the CMB reconstruction. This will be better discussed 
in the next section, where we will focus our analysis on the CMB 
component. Variations are very small also for the reconstructed 
thermal dust and free-free components. This is due to the fixed fre- 
quency scaling and to the use of prior maps. 

As expected, we have a substantial spread in the reconstructed 
synchrotron component, reflecting the wide range of synchrotron 
spectral indices used for the separation. This component is found to 
be generally sub-dominant compared to the free-fr ee, except, pos- 
sibly, in the K band, in agreement with the results bv lHinshaw et al.l 
(2006). 

The results on the anomalous emission components are also 
quite stable. If we force this component to have a power-law spec- 
trum (model M3) we find that it dominates over both synchrotron 
and free-free in all the WMAP channels. The anomalous emission 
for M4 is found to dominate in the Ka, Q and V bands. Models 
M2 and M5 are consistent with each other, as the summed inten- 
sity of thermal dust and anomalous emission of model M2 is close 
to that of the "total dust" emission of model M5. Moreover, they 
yield similar intensities for the remaining components. As previ- 
ously stressed, models M2 and M 5 yield results on the anomalous 
emission which are consistent with lDavies et al. 

5.1 Quality tests 

Before proceeding with the analysis of the results we try to evalu- 
ate the goodness of the decomposition for each model. A standard 
method is to analyze the residual between the data and the recon- 
structed sources combined by means of the estimated mixing ma- 
trix: 

r = X - Hs. (22) 

The analysis of the residuals allows to check if, given a cer- 
tain mixing matrix, our algorithm succeeds in finding a plausible 
decomposition of the data. In the case of a perfect reconstruction of 
the channels, r contains maps of pure noise, characterized by Gaus- 
sian statistics. Foreground residuals bring in non-Gaussianities. A 
visual inspection of the residuals show that some non-Gaussian fea- 
tures are present in the maps. In particular, we can discern traces of 
anisotropic noise, point sources and diffuse emissions, particularly 
at low Galactic latitudes. The first two features are related to effects 
that are not accounted for in our analysis. The latter is more inter- 
esting, because it also reflects our modeling and estimation errors. 

Each model shows the worst contamination in the K band, 
even if models including the anomalous emission are generally 
more clean (see for example the comparison between Ml and M2 
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Figure 8. Frequency dependence of the rms fluctuations, in thermodynamic temperature, of the components reconstructed with each model. The shaded ai'eas, 
visible only for components whose spectral parameters are allowed to vary, show the dispersion obtained from the 30 reconstructions per model, and the lines 
shows the mean spectra. 



in Fig.|9](. The "residual excess" at low latitudes then decreases at 
higher frequencies. The models Ml and M4 are generally worse; 
however this strongly depends on the synchrotron spectral index 
used for each separation. 

More quantitative analyses of the residuals (power spectra, 
standardized moments of the distributions, etc.) mainly confirm 
these general findings, but turned out to be not very effective in 
discriminating among the models, as the results are comparable. 

This means that, in our case, each model allows a satisfactory 
decomposition of the data, at least for a subset of synchrotron spec- 
tral indices. Neverthless, some of the models could still give uncor- 
rect descriptions of the data. The analysis of the residuals, infact, 
is not informative on the reliability of the individual reconstructed 
components. For example, they may be affected by aliasing, which 
does not show up in the residuals. To investigate this possibility, 
we applied other quality tests directly based on the reconstructed 
components. 

5.1.1 Correlation among the components 

Independent observations have highlighted positive correlations of 
varying strength among synchrotron, free-free and dust emissions, 
ensuing from the physical processes producing them. Therefore, 
the reliability of the recovered components can be tested by inves- 
tigating their mutual correlations. We have worked with pixels of 



area 0.84 deg^ (Nside = 64). We trimmed all pixels pre-processed 
through the apodized mask of Fig. [7] as well pixels within a ra- 
dius of 1° around each point source listed in the WMAP three-year 
catalogue. We also excluded a region with a radius around 15° 
around the Gum Nebula, as in this region the separations are al- 
ways very noisy. High Galactic latitude regions (|ti| > 50°), where 
foregrounds are weak and therefore recovered with low S/N ratios 
have also been excluded from the analysis. Around 3 • 10* pixels 
then remain. 

For each model we have computed the correlation coefficients, 
r, among the recovered Galactic components for each of thirty sep- 
arations. Whenever one of those coefficients was found to be nega- 
tive, the corresponding separation was discarded. The percentages 
of discarded separations for each model are reported in Table [3] 
They are generally low 10%), except in the case of M3 (30%). 

All the models including the anomalous component feature 
very strong positive correlations of it with thermal dust (r ~ 0.8), 
somewhat less strong but still highly significant positive correla- 
tions with synchrotron (r ~ 0.6), and weak positive correlations 
with free-free (r ~ 0.25). 

For the large pixel number, A'^ ~ 3 x 10*, we are dealing with 
the correlation coefficient, r, for samples of uncorrelated quantities 
has a normal distribution with zero mean and Or = l/y/W. CMB 
and Galactic foregrounds are intrinsically uncorrelated; therefore, 
the presence of a statistically significant correlation among them 



WMAP Syr data with the CCA 9 



Q bond 







j 










j 








































]- 




















.-■ a i 


f 1 


■'•i- 






?* 




'■Kl ■ ■ 1 

At* ]>. J-V ^\ 






jgi i : 














■ 1 


-7.20 










-7.30 


-dl ■ 








-Dl 



0.100 



0.010 



0.00 




2.80 




V band 



Figure 9. Residual at K band for model Ml (left) and M2 (righ) assuming 
jSs ~ 3. in two regions of ~ 40X40 degrees. Central coordinates of the 
patches are Z = 0°, 6 = 15° (up) and I = 140°, b = 15° (down). Units 
are mK thermodynamic. 



Table 3. Percentage of discarded separations because of negative correla- 
tion coefficients. 
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Figure 10. Absolute values of the correlation coefficient between the recov- 
ered CMB and the total foreground map for the Q (up) and V (down) bands 
as a function of the synchrotron spectral index. Values above dashed line 
correspond to congelations formally significant at Scr. 



model 


discarded separations 


Ml 


0% 


M2 


6% 


M3 


30% 


M4 


10% 


MS 


3% 



indicates a poor component separation. Thus, a good criterion to 
evaluate the quality of the separations is to compute the correlation 
coefficient between the reconstructed CMB and foreground maps. 
The latter have been obtained as the sum of the reconstructed com- 
ponents scaled to the Q and V bands. As in these bands the fore- 
ground contamination has a minimum, those correlations are ex- 
pected to be low and, hopefully, not statistically significant. Once 
computed for each separation and each model, the correlation coef- 
ficients can be used as figures of merit to compare different models 
and, for a given model, to identify the best separations. 

In Figure[TO]we show the absolute values of the correlation co- 
efficients as a function of the assumed synchrotron spectral index. 
The horizontal dashed line is the 3(t significance level. 

Model M3 always features highly significant correlations (r > 
0.1) between the CMB and the reconstructed foreground contribu- 
tion in both bands, much stronger than found for the other models. 
This result, together with the high number of samples we had to 
discard due to aliasing effects, indicates that this model is not able 



to correctly describe the data. It was therefore discarded together 
with model M4 which also yields values of jrj always well above 
the 3(T limit, although somewhat below those found for M3. 

In the case of models Ml, M2 and M5, the test shows that the 
correlations are not statistically significant for well defined ranges 
of the synchrotron spectral index. Model Ml requires relatively low 
values of j3a, below that expected, in the WMAP frequency range, 
from the locally measured energy spectrum of relativistic electrons 
which would yield ~ 3 (Bandav & Wolfendale 1990, 1991), as 
indeed found for models M2 and MS. 

We have built reference maps of each reconstructed compo- 
nent for each of the 3 surviving models by averaging those recon- 
structions with foreground-CMB correlations at less than Scr level. 
For model M2 we find somewhat different allowed ranges of (is 
for the Q and the V bands, and we have kept the reconstructions in 
the union of the 2 ranges, i.e. with 2.9 ^ /3s ^ 3.1. In Table |4]we 
list the correlation coefficients of the reference CMB map for each 
model with the foreground templates. The correlations computed 
with the minimal mask are formally statistically significant, indi- 
cating the presence of residual foreground contamination. The ma- 
jor contributions to the correlation coefficients come from regions 
close to the Galactic plane. If we adopt the WMAP kp2 mask, in 
place of our minimal mask, the correlation coefficients decrease, 
sometimes going close to or below the 3cr significance level, es- 
pecially for free-free. In Table |5] we show the correlations of the 
reconstructed synchrotron, dust and free-free components with the 
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Table 6. Correlation coefficients of the anomalous emission (M2) and of 
the total dust component (MS) with all the foreground templates. 
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M2 
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Figure 11. CMB map reconstructed with model M2. 



Table 4. Correlation coefficients of the reconstructed CMB components 
with the foreground templates: minimal mask, (kp2 mask) 



template 


Ml 


M2 


MS 


thermal dust 


0.027 (0.015) 


0.029 (0.023) 


0.025 (-0.023) 


synchrotron 


-0.054 (-0.046) 


-0.016 (-0.012) 


-0.015 (-0.010) 


free-free 


-0.035 (-0.021) 


-0.035 (-0.0003) 


-0.044 (-0.011) 




Multipole moment 

Figure 14. Full-sky power spectra of the foreground templates at low mul- 
tipoles (arbitrary normahzations) 



corresponding templates. The match between the prior and recon- 
structed maps is good but not perfect, as expected given the spa- 
tial variations of the spectral properties. Finally, Table |6] gives the 
correlation coefficients between the anomalous emission recovered 
with model M2, and the thermal dust plus anomalous emission re- 
covered with model M5 (which treats them as a single component 
with a complex spectrum), and the foreground templates. Interest- 
ingly, in the case of model M2, which does not impose a priori any 
correlation between thermal dust and anomalous emission, we find 
a tight, but not perfect, correlation among the 2 components. The 
anomalous emission also correlates strongly with synchrotron, and 
more weakly with free-free. 

The final map of the anomalous component recovered with 
model M2 is shown in Fig. [12] where the dust and synchrotron 
templates are also shown for comparison. Traces of imperfect com- 
ponent separation can be discerned by eye: negative imprints in 
correspondence of the Gum Nebula, of the Orion region and of 
the North Polar Spur, as well as several point-like sources. Nev- 
ertheless, this may be the first, albeit preliminary all-sky map of 
this component. The similarity with the thermal dust map is evi- 
dent; we recall that no prior correlation of the anomalous emission 
with thermal dust was imposed. Differences with the spatial dis- 
tribution of synchrotron are also clear: the anomalous emission is 
relatively stronger towards the outskirts of the Galaxy, where the 

Table S. Correlation coefficients of the synchrotron, thermal dust and free- 
free maps with the con'esponding templates. We recall that for MS the dust 
component includes both the thermal and the anomalous contributions. 



template 


Ml 


M2 


MS 


thermal dust 


0.966 


0.967 


0.883 


synchrotron 


0.865 


0.858 


0.851 


free-free 


0.957 


0.984 


0.984 



first application of CCA found flat "synchrotron" spectral indices, 
and less prominent near the Galactic centre. In Fig.[T3]we show a 
detail of the comparison between the thermal dust and anomalous 
emi ssion maps recovered with M2. The area we show is the Region 
6 of lOavies et al.l ^20061) . chosen to be dust dominated and poorly 
contaminated by free-free and synchrotron. In agreement with their 
results, the anomalous emission is detected in this region, and ap- 
pears to be dust-correlated. 



6 THE CMB POWER SPECTRUM 

The next question we wish to address is: to what extent do the dif- 
ferent foreground models affect the estimates of the CMB power 
spectrum on large scales, where possible inconsistencies with pre- 
dictions of the standard cosmological paradigm have been re- 
ported? 

We will focus our analysis on I ^ 70, where the effect of 
fluctuations due to unresolved point sources, non included in our 
study, can be neglected. On large scales, however, the CMB maps 
reconstructed through matrix inversion of the WMAP data comple- 
mented with the three template foreground maps, keep trace of the 
spurious structure due to latter maps being mosaics of observations. 
This is clearly visible in Fig. [14] showing the power spectra of the 
templates at low multipoles. While for Z ^ 10 the power spectra are 
smooth, at lower multipoles there are conspicuous spikes, particu- 
larly for the free-free and the synchrotron. To avoid this problem, 
the CMB power spectrum for / < 10 was computed using only 
the WMAP data for the Wiener Filter reconstruction of the CMB, 
still using the same ranges of optimal spectral parameters for each 
model. For larger Vs. the inversion is carried out including the fore- 
ground templates. T he error analysis, carried out as described in 
iBonaldi et al.l ( l2006h . confirms that this approach significantly de- 
creases the errors on the power spectrum at low Ts. The power 
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Figure 12. Map of the anomalous emission recovered witli model M2 in the K band compared with the dust (850 GHz) and the synchrotron (408 MHz) 
templates. Minimum and maximum value of the colour scales are given by the mean of the map ± its standard deviation. 



spectra obtained with the two methods are very close to each other 
for I ~ 10, so that they smoothly join and the choice of the bound- 
ary multipole needs not to be fine tuned. 

In these calculations we used two sky masks. The first one 
(minimal mask) excludes all pixels having non unitary value in the 
apodized mask used to preprocess the data (Figure|7]l; the cut region 
amounts to ~ 10% of the sky. The second is the previous one mul- 
tiplied by the WMAP kp2 mask; it excludes ~ 20% of the sky. The 
latter is less liable to residual foreground contamination, while the 
former is less liable to biases due to incomplete sky coverage. The 



power spectrum was binned with the same scheme ad opted by the 
WMAP team and we applied the MASTER approach jHivon et al.l 
2002) to get an unbiased estimate. 

The results for each model are shown in Fig. [15] where we 
have plotted the power spectra computed with the minimal mask 
for Z ^ 8, those computed with the other mask for higher mul- 
tipoles. Again the results are weakly dependent on the choice of 
the mask, but the minimal mask yields smaller oscillations at low 
Vs. We find generally a good agreement with the WMAP power 
spectrum, shown by the solid line. There are however significant 
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differences among the various models at tiie lowest multipoles, es- 
pecially evident for I = 2. Regardless to which are the actual mod- 
els we are exploiting, this result shows that the hypotheses on the 
Galactic components made to perform component separation can 
be painful for the CMB power spectrum at larger scales. Thus, big 
error bars are required to take into account for possible biases at 
these multipoles. As an estimate of uncertainties associated to fore- 
ground modeling, we can take the spread of our three CMB power 
spectra. 

In the upper panel of Fig.[T6]we compare our results with the 
WMAP three year power spectru m and with the best fit ACDM 
model based on WMAP data only jSpergel et alj|200^ . The mean 
quadrupole moment coming out from our analysis is higher than 
the WMAP estimate, and the difference with the prediction of the 
ACDM model is within the uncertainty due to cosmic variance. 
The spread of estimates of the quadrupole moment from different 
models is of about ±200/iK'^. 

Summing in quadrature the cosmic variance and the modeling 
errors, we find no large scale power spectrum "anomalies" signifi- 
cant at ^ 1.5(T, except for the excess power at I ~ 40, which is sig- 
nificant at ~ 4(T. On the other hand, as show by Fig.[T7] the North- 
South asymmetry at Z ^ 20 stands, independently of the adopted 
foreground model. Those power spectra, computed for each model 
for the North and for the South hemisphere separately, were ob- 
tained as previously described, except for the different binning 
scheme adopted. 



7 CONCLUSIONS 

We have performed the component separation on WMAP maps ex- 
ploiting CCA technique to estimate the mixing matrix directly from 
the WMAP data, t aking also into account the complementary in- 
formation from the lHaslam et al.l jl982h 408 GHz map, from the 
extinction corrected Hq map jPickinson et alj2003l), and from th e 
combined COBE/IRAS thermal dust maps ( lFinkbeineretal.ll999h . 
The Ha map is exploited also to subtract the free-free contribution, 
which is subs t antial in some regions of the Galactic plane, from the 
iHaslam et al. I( ll982l) map, thus obtaining a pure synchrotron tem- 
plate. 

The CCA is particularly useful to deal with possible additional 
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Figure 15. Binned CMB power spectra obtained from all models compared 
to the WMAP thi'ee-year power spectrum 
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Figure 16. Upper panel: best fit ACDM model(smooth solid line) com- 
pared to the WMAP CMB power spectrum (dashed line) and to the mean 
final power spectrum resulting from our analysis (open diamonds; the er- 
ror bars show the spread of results for different models). The shaded area 
shows the cosmic variance. Lower panel: uncertainties associated to fore- 
ground modeling (dotted line) compai'ed to the cosmic variance (solid line). 
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Figure 17. North-south asymmetry. The bins adopted have central multi- 
poles / = {3, 8, 13, 18, 23, 28, 33, 38} and equal width A/ = 5 



foreground components, for which little or no prior information ex- 
ists. An application assuming that only the canonical Galactic emis- 
sions (synchrotron, free-free, and thermal dust) are present high- 
lights the widespread presence of a spectrally flat "synchrotron" 
component (Fig.[Tll, largely uncorrelated with the synchrotron tem- 
plate (compare Fig. |2] with the synchrotron template in Fig. I12t . 
suggesting that an additional component is required. 

We have then tested various spectral shapes for such compo- 
nent, namely a power law (model M3) as expected if it is a flat syn- 
chrotron emission due to enhanced magnetic fields in dusty star- 
forming regions, or a parabola in the log S-log u plane (model 
M2), or a grey body, as proposed by Tegmark et al. (2000; model 
M4). In all these cases, the spatial distribution of this component 
was left totally unconstrained, but we have also tested the possibil- 
ity that it is distributed exactly as thermal dust (model MS). 

The CMB maps yielded by models M3 and M4 turned out to 
be substantially correlated with foreground maps, indicating a un- 
acceptably large residual contamination. These models were there- 
fore discarded. Conversely, in the case of models M2 and MS the 
correlations were found to be not statistically significant for values 
of the synchrotron spectral index, f3s — 3, close to the expectations 
from the slope of the locally measured energy spectrum of rela- 
tivistic electrons (Banday & Wolfendale 1990, 1991). The data do 
not allow us to clearly discriminate among these models; however, 
some interesting indications emerge: 

• The additional component turns out to be always tightly cor- 
related with thermal dust although, if its spatial distribution is un- 
constrained, the correlation is not perfect (co rrelation coefficien t 
r ~ 0.8). This reminds us the suggestion bv Davies et al. I j2006l) 
that, if this component is due to spinning dust, it should be better 
correlated with the small grains dominating the mid-IR emission 
than with the big grains dominating at far-IR to sub-mm wave- 
lengths. The map of this component yielded by model M2 could 
then constitute its first, albeit preliminary all sky map. 

• The additional component is well correlated with synchrotron 
(correlation coefficient r ~ 0.6) and more weakly correlated with 
the free-free, although the free-free is strongly correlated with ther- 
mal dust (correlation coefficient r — 0.4 — 0.6). The synchrotron 
is also well correlated with thermal dust (correlation coefficient 
r = 0.7 — 0.8), and more weakly with the free-free (correlation 
coefficient r ~ 0.2). 

• If only the standard Galactic emissions (synchrotron, free-free 



and thermal dust) are taken into account, we obtain acceptable sep- 
arations adopting a flat synchrotron spectral index (f3s ~ 2.8). 
However, this value of Ps is not supported by the CCA analysis 
and is flatter than expected in the WMAP frequency range (Banday 
& Wolfendale 1990, 1991). 

• The CMB maps we obtain have low foreground contamina- 
tion even at low Galactic latitude, as demonstrated by the correla- 
tion with the foreground templates. The inclusion of the anomalous 
emission decreases the synchrotron contamination. 

We have used the MASTER approach (Hivon et al. 2002) 
to obtain unbiased estimates of the CMB power spectrum up to 
I — 70 from maps reconstructed on the basis of each foreground 
model. The results are generally consistent with the WMAP ones, 
although we find a substantial spread of values for the quadrupole 
moment. For model M2 the quadrupole moment is fully consistent 
with the expectation from the "concordance" cosmological model. 
The spread of power spectrum estimates gives us a measure of un- 
certainties associated to foreground modeling. Such uncertainties 
are smaller than the cosmic variance, but nevertheless significant, 
for I ^ 30, and approach it for larger Vs. Combining the cosmic 
variance with the modeling errors, we find that the quadrupole 
amplitude is less than la below that expected from the standard 
ACDM model. Also the other reported deviations from model pre- 
dictions are found not to be statistically significant, except for the 
excess power at I ~ 40. On the other hand, the North-South asym- 
metry is found to be unrelated to foreground models. The derived 
CMB power spectra are remarkably stable: even those yielded by 
the discarded models are within the range of Fig. 1 161 
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